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ABSTRACT 

Using deep Subaru/FOCAS spectra of 34 Hn regions in both the inner and outer 
parts of the extended ultraviolet (XUV) disc galaxy NGC 4625 we have measured an 
abundance gradient out to almost 2.5 times the optical isophotal radius. We applied 
several strong line abundance calibrations to determine the H n region abundances, 
including R 2 3, [Nil] /[On], [Nil] /Ha as well as the [Om] A4363 auroral line, which we 
detected in three of the Hn regions. We find that at the transition between the inner 
and outer disc the abundance gradient becomes flatter. In addition, there appears to 
be an abundance discontinuity in proximity of this transition. Several of our target 
Hn regions appear to deviate from the ionisation sequence defined in the [Nil] /Ha vs. 
[Oin]/H/3 diagnostic diagram by bright extragalactic Hn regions. Using theoretical 
models we conclude that the most likely explanations for these deviations are either 
related to the time evolution of the H n regions, or stochastic variations in the ionising 
stellar populations of these low mass H n regions, although we are unable to distinguish 
between these two effects. Such effects can also impact on the reliability of the strong 
line abundance determinations. 

Key words: galaxies: abundances - galaxies: structure - galaxies: ISM - galaxies: 
individual: NGC 4625. 



1 INTRODUCTION 

Extragalactic Hn regions act like galactic buoys, highlight- 
ing sites of recent or ongoing star formation. It has long 
been known that H II regions are not just found in the ac- 
tive star forming discs of galaxies, but also at large galac- 
tocentric radii, as shown, for example, by [Ferguson et al.| 
( 1998b I and Lelievre & Roy ( 2000 1. These observations have 



been bolstered by the discovery of young B-type stars at 
large galactocentric radii (Cuilla ndre et aL]|2001| |Davidge] 
20071. However, the idea of star formation occurring be- 



yond the optical edge of a galaxy (commonly defined by the 
25th magnitude B-band isophotal radius, R25) was revolu- 
tionised with the advent of the Galaxy Evolution Explorer 
(GALEX) satellite ( |Martin fc GALEX Science Team|2003 ). 
The use of GALEX as part of the Nearby Galaxies Survey 
(Bian chl et al.|[2003 ') revealed numerous galaxies where the 
UV emission extends well beyond the optical edge. These 
galaxies have since been termed 'extended ultraviolet disc' 



(XUV disc) galaxies (Thilker et al. 
et al. 12005) . 
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The GALEX satellite, in conjunction with optical and 
infrared observations, has given us some major insights into 
the nature of these XUV discs. We now know that these 
structures are rather common, with an incident rate > 30% 
in spiral galaxies (Zaritsky & Christlcin 2007 Thilker et al. 
|2007[). XUV disc s are complex and varied systems, divided 



by Thilker et al. ( 2007 1 into two types. Type 1 XUV discs are 



highly structured, displaying UV-bright, yet optically dim, 
knots. Type 2 XUV discs, which are predominantly found 
in late- type spirals, show an extended UV-bright, but once 
again optically faint, region beyond the galaxy edge. Many of 
the UV knots correspond to structures in the H I distribution 
( Gil de Paz et al.|2005| |Thilker et al.|2005[), and some have 



Hn region counterparts (Gil de Paz et al. 2005 Goddard 



|et al.||2010"l ). The UV and Ha fluxes of these Hn regions 
indicate a mass of the underlying stellar population between 
10 3 and 10 4 M®. These studies also show that ionisation is 
dominated by only a handful or a single Q-type star (|Gil de| 



Paz et al.|2005)|Bresolin et al.|2009a||Goddard et al.|2010 | 

The ages of the UV complexes in XUV discs appear to range 
between a few Myr and several hundred Myr ( jZaritsky fc 
Christlein|2007| |Dong et al.|2008 l. 
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Even though the mean surface gas density at extreme 
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galactocentric radii falls below the critical threshold for star 



formation (Kennicutt 1989 1, local density perturbations may 
still promote star forming activity in the outer discs of spiral 
galaxies JMartin & Kennicutt 2001). This is supported by 



Table 1. Hn region sample. 



models of the propagation of spiral density waves into the ex- 
tended discs ( Bush et al.|2008 l. However, other mechanisms 



may also play a part. Gil de Paz et al. (I2005J) , for exam- 
ple, suggested that an interaction between NGC 4625 and 
NGC 4618 acted as a trigger for star formation at extreme 
radii in NGC 4625. lElmegreen & Hunterl (120061) proposed 



a range of other mechanisms, including stellar compression 
and turbulence compression among others. 

Hn regions in XUV discs are ideal tracers of star- 
formation at extreme radii, but they also encode information 
about the kinematics of these structures ( |Christlein fc ZariFj 
sky||2008 1 and the chemical evolution of these extreme star 



formation environments ( |Gil de Paz et aL 2007 Bresolin 



|et al.|[2~009a| ). Previous studies of the abundances of XUV 
disc H II regions have insofar been rather limited. |Ferguson] 
eTal] (|1998a|) studied the galaxies NGC 628, NGC 1058 and 



NGC 6946, but their sample was limited to nine Hn regions 
beyond R25 . More in depth studies of H 11 region abundances 
in XUV discs have focused on two galaxies in particular. |Gil| 
de Paz et al. (2007 hereafter referred to as G07) studied 



31 Hll regions in the XUV discs of M83 and NGC 4625, 
finding in both cases an exponential abundance gradient ex- 
tending into the outer disc. However, a deeper study of M83 
by Bresolin et al. (2009a I, containing 49 Hn regions out to 



2.6 R25, found a marked discontinuity in the abundance gra- 
dient at the optical edge of the galaxy, together with a flat 
gradient in the XUV disc. 

NGC 4625 and M83 have been the most studied XUV 
disc galaxies, partly due to their close proximity. These 
galaxies possess some of the most XUV discs known, ex- 
tending out to 4 R25 in both cases ( |Thilker et al.|2005| [Gn] 
de Paz et al.|[2005 1. They share similar trends in their Ha 
and UV radial profiles, with a turnover in the Ha profile 
near R25, yet a smooth far-UV profile extending into the 
XUV disc ( |Goddard et al.||2010[ ). However, M83 displays a 
large number of UV knots in the XUV disc in association 
with structured Hi emission, whilst the lower luminosity 
NGC 4625 has a more diffuse UV emission accompanying 
a low surface brightness optical component (Thilker et al 



2007) 



In this paper we present a detailed chemical abundance 
study of the Hn regions in NGC 4625, including targets 
across the XUV disc not previously observed by G07. We 
have sampled Hn regions in both the inner and outer disc 
out to ~ 2.5 R25, and study their oxygen abundances based 
on several strong nebular emission line indicators. We first 
describe our observations and the data reduction in Ej2] In 
Sj3]we discuss the oxygen abundances that we obtained and 
compare our results to those of previous datasets. ij4] dis- 
cusses the nature of Hn regions in the XUV disc, consid- 
ering the effects of an ageing population and stochasticity 
on the observed nebular line ratios. We briefly discuss our 
findings in the context of extragalactic star formation in fj5| 
before we finally summarise our conclusions. Throughout 
this paper we assume a distance to NGC 4625 of 9.5 Mpc 
( Kennicutt et alT||2003 1, an inclination angle i = 27° and a 
position angle of the major axis 9 = 1 50° (|Bush fc Wilcots 
2004 ). The R25 value of 66 arcsec, from de Vaucouleurs et al. 
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(1) Identification number; a compound number (e.g. 131.228) in- 
dicates an Hn region observed with both our multi object masks; 
ID numbers with a suffix (e.g. 116c) indicate that more than one 
H II region was extracted from the same slit. XUV numbers identify 
the objects studied by G07. (2, 3) Right Ascension and Declina- 
tion (J2000). (4) Galactocentric distance in units of the isophotal 
radius R 25 (R25 = 66" = 3.04 kpc). 



( 1991 1, corresponds to 3.04 kpc. We adopt a solar metallicity 
value 12 + log(O/H) =8.69 ( |Asplund et al.|2009 l. 



2 OBSERVATIONS AND DATA REDUCTION 

Our H 11 region selection was based on deep Ha images taken 
using the Faint Object Camera And Spectrograph (FOCAS, 
Kashikawa ct al. 2002| at the Subaru telescope situated on 
Mauna Kea in Hawaii. The field of view of the instrument 
(6 arcmin) encompasses the entire XUV disc of NGC 4625. 
A continuum-subtracted Ha image was constructed using a 
900 s on-band image in combination with a 120 s continuum 
image. We designed two multi-object slit masks, attempting 
to cover as wide a range in galactocentric radius as possi- 
ble, and including several Hn regions which were part of 
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the spectroscopic study by G07, to enable a comparison of 
identical Hn regions. 

The spectra were obtained with FOCAS on the night 
of March 22, 2009, using 1" (= 46 pc) slits in combina- 
tion with the 300R grism (used in second order) to cover 
the blue part of the spectrum (approximately from 3400 A 
to 5300 A) and the VPH650 grism for the red part of the 
spectrum (5300 A to 7700 A). The spectral resolution, esti- 
mated from the gaussian FWHM of lines in the arc lamp 
frames, is approximately 5.5 A. We exposed for a total time 
of 2 hours in the blue and 80 minutes in the red for each 
mask. Seeing conditions varied through the night from a 
poor 2" to a more favourable 0'.'8 on occasion. The varia- 
tions in the atmospheric conditions occurred mostly during 
the observations taken with the first mask, while the seeing 
stabilized around 1" for the remaining part of the night. To 
flux-calibrate our spectra we took long slit spectra of four 
standard stars, GD71 and G191-B2B at evening twilight, 
and HZ44 and BD+33d2642 at the end of the night. The 
data were reduced using standard iraiQ routines. 

We obtained spectra for 34 Hn regions from the two 
slit masks combined. We list their positions and galacto- 
centric distances in Table [l] Our targets are marked on a 
continuum-subtracted Ha image and a GALEX ultraviolet 
colour composite image of NGC 4625 in Fig. [I] Four of the 
Hn regions were observed in both masks. A few of the slits, 
identified by a letter following the slit identification number 
in Table [I] included multiple targets. 

We measured the spectral line strengths by direct inte- 
gration of the flux measured between two continuum regions 
interactively selected on each side of the line. The reddening 
correction was done using the extinction curve published by 
|Qsterbrock| ( |198"9| , and assuming an intrinsic H7/H/? ratio 
of 0.47, which holds at T e = 10 4 K. We were unable to use 
the Ha/H/3 ratio because these two Balmer lines appear sep- 
arately in the red and blue spectra, respectively. Therefore, 
the theoretical ratio Ha/H/3 = 2.86 was adopted to renor- 
malise the line fluxes in the red spectral range, after we had 
made the reddening correction. The line flux ratios that in- 
volve emission features observed separately in the blue and 
red spectra could potentially be affected by the variations 
in the seeing conditions mentioned earlier. Fortunately, the 
spectra for our second mask were obtained with an approxi- 
mately constant image quality (~ 1"). Moreover, since a few 
of the targets were included in both masks, we could verify 
that the flux ratios did not suffer from changes in the seeing. 

We list the reddening-corrected line fluxes of our sam- 
ple, in units of H/3 = 100, in Table [2] The quoted errors 
account for measurement, flux calibration and flat fielding 
errors. We point out that the H/3 flux (column 7) underesti- 
mates the total line flux (hence the ionizing luminosity, see 
j ]4.2.1[ |, due to the unavoidable slit losses. 

We have characterised our Hn region sample using the 
familiar BPT (Bal dwin et al.|1 981) diagnostic diagrams, by 
plotting in Fig. |2| [ N n]/Hq against [Om]/H^ an d [Sii]/Hq 
against [O ill] /H/3 ( |Veilleux fc Osterbrock||1987| >. We repre- 



sent our sample of H 11 regions with orange triangle symbols, 
and we include comparison data of mostly bright extragalac- 
tic H 11 regions drawn from the literature: 

1. 99 Hn regions measured in 20 spiral and irregular 



galaxies by McCall et al. (1985 dark grey dots). 
2. 198 H 11 regions in 13 spiral galaxies from | van Zee et aL] 



|1998[ grey dots). 
3. 17 Hn regions from three late-type spirals (NGC 628, 



NGC 1058 and NGC 6946) studied by |Ferguson etaL] p998a 
light grey dots). 
4- 41 Hn regions in the spiral galaxy M101 observed by 



Kennicutt & Garnett (1996 black dots). 

5. 16 H 11 regions from the inner disc of M83 from |Bresolin| 
& Kennicutt (2002) dark blue squares). 

6. 49 H II regions from the M83 sample of Brcsolin et al. 
( 2009a] light blue squares), mostly composed of outer disc 
Hn regions. Note that the two objects which appear above 
the theoretical boundary in Fig. [2] were identified as super- 
nova remnants. 

Hn regions in the inner disc of NGC 4625 are 
further identified by green dots, while the outer disc 
Hn regions are marked either with cyan dots, in- 
dicating [0 m] AA4959,5007/H,fl > 0.25, or red dots 
([Oiii] AA4959, 5007/H/3 < 0.25). This arbitrary distinction 
helps us to distinguish these two subsets. In fact, in the 
[Nil] /Ha vs. [Oiii] /H/3 diagram not all of the inner disc Hn 
regions sit on the tight sequence defined by bright Hn re- 
gions: those marked with red dots appear to lie below and/or 
to the left of the main ionisation sequence. On the other 
hand, the [Sll]/Ha versus [O ill] /H/3 plot shows a similar 
ionisation sequence for objects in the inner and extended 
disc of NGC 4625, which is consistent with the ionisation 
sequence defined by bright extragalactic Hn regions. These 
results will be further discussed in S)4] 

Eight of our H 11 regions were also studied by G07 (iden- 
tified in Table [I] by the xuv numbers used in their paper) . 
The H/3 fluxes measured by G07 are systematically higher 
than ours by a mean factor of 1.45, as expected from their 
50% larger slit widths (l'/5 vs. l'.'O), and thus it appears 
that the poor seeing conditions that we experienced did not 
significantly affect our flux measurements. Our [O ill] /H/3 
ratios for 6 H 11 regions in common (not all of the G07 had 
a measured [O ill] line) are systematically 22% lower, on 
average, for which we do not have an explanation. Of the 
only four targets for which the intensity of the [0 11] A3727 
line could be compared, we found good agreement in two 
cases (our regions 118 and 122), a 34% lower flux ratio in 
a third (region 103), which is mostly due to our lower as- 
sumed reddening correction, while in the case of region 130 
our [0 11] /H/3 value is almost a factor of 2 lower. In this last 
case, the reddening correction cannot explain the large dis- 



1 IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under co-operative agreement with 
the National Science Foundation. 



crepancy. It is well-known (e.g. Bresolin et al. [2005) that 
the [0 11] A3727 line is susceptible to systematic errors in 
nebular optical spectra, as a result of the increased uncer- 
tainty in the flux calibration of the blue end of the spectrum, 
smaller throughput at these wavelengths, and errors in the 
reddening corrections. We conclude by noting that most of 
the spectra shown in Fig. 4 of G07 are very noisy in the 
region around 4000 A, and in several cases even the H/3 line 
seems to be detected with a poor signal-to-noise ratio (un- 
fortunately G07 did not publish error estimates for their line 
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Figure 1. Left: Continuum-subtracted Ha image of NGC 4625, obtained from our Subaru FOCAS imaging. Right: GALEX UV colour 
composite image of NGC 4625. Our Hll region sample is shown by the orange circles marked with the object IDs. 




Figure 2. BPT diagnostic plots. Left: [Om]/Bfi vs. [Nil] /Ha. Right: [OiiiJ/H/3 vs. [Sn]/Ha. Our sample of Hll regions in NGC 4625 
(triangles) has been plotted together with published Hll region datasets. Objects in the inner disc are further identified by green dots, 
while the outer disc H II regions are marked either with cyan dots, indicating [O lll]/H/3 > 0.25, or red dots for [O IIl]/H/3 < 0.25. The blue 
squares show the Hll regions observed in the other prototypical XUV galaxy, M83, with light blue squares representing the outer disc 
sample of [Bresolin et aL] | |2009a| , and dark squares the inner disc Hll region s of|Bresolin &; Kennicutt] (|2002fr. The do ts show previous 
studies of bright inner disc H II regions, including the samples of |McCall et al.| | |1985| dark grey dots) and|van Zee et al.|(|1998| grey dots). 
The light grey dots represent the H II regions in NGC 628, NGC 1058 and NGC 6946 studied by |Ferguson et al.| | |1998a| l' The M101 
sample of |Kennicutt fc Garnett|( 1996|l is represented by the black dots. The red curves represent the upper boundaries for photoionised 
nebulae defined by Kcwlcy et al. (2006J. The dashed lines in the left panel serve to schematically isolate objects that appear to deviate 
from the main ionisation sequence, towards the lower left of the diagram. 
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Table 2. Reddening-corrected lines fluxes. 
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30 ± 1 


105 ±4 


74 ±2 


5 ±0.3 


18.9 


0.00 ±0.13 


-0.28 


-0.43 


214 


149 ±9 


28 ±2 


96 ±4 


83 ±3 


5 ±0.9 


14.4 


0.00 ±0.21 


-0.19 


-0.47 


215 


225 ± 17 


54 ±2 


112 ± 2 


149 ±2 


8 ±0.2 


5.7 


0.40 ±0.14 


-0.30 


-0.41 


216 


230 ± 32 


256 ± 13 


68 ±3 


104 ±4 


10 ±2.1 


2.6 


0.00 ±0.19 


-0.53 


-0.62 


217 


192 ±8 


76 ±3 


73 ±3 


58 ± 2 


8 ±0.5 


31.7 


0.00 ±0.13 


-0.42 


-0.59 


221a 


359 ± 98 


164 ± 19 


43 ±2 


75 ±3 




1.2 


0.27 ± 0.44 


-0.92 


-0.82 


221b 


353 ± 43 


75 ±6 


31 ±4 


84 ±6 




0.9 


0.00 ±0.32 


-1.05 


-0.96 


223 


357 ±88 


55 ±6 


39 ±2 


87 ±3 




2.9 


0.58 ±0.39 


-0.95 


-0.86 



Line fluxes are normalized to H/3 = 100. The H/3 flux in col. 7 is uncorrected for potential slit losses. 



ratios, so this statement is based on a visual assessment of 
their data). Thanks to longer exposure times (2 hours vs. 
40 minutes) with a larger telescope aperture (Subaru 8 m 
vs. Palomar 5 m), our spectra are, as should be expected, of 
superior quality. We are therefore confident in the reliabil- 
ity of our emission line data and the estimated uncertainties 
presented in Table [2] 

3 NEBULAR ABUNDANCES 
3.1 Direct Abundances 

Ideally, the derivation of the chemical abundances of ionised 
nebulae requires first a measurement of the electron tem- 
perature T e , since the line emissivities strongly depend on 
it. The classical method to directly obtain T e involves the 
simultaneous measurement of the auroral [O in] A4363 line 
and of the nebular [O in] AA4959, 5007 lines. This proves to 
represent a challenge with the increase of the metallicity, 
because the cooling mechanisms become highly efficient and 



[O in] A4363 too weak to observe. However, having a few 
A4363 detections is important, since they allow the deriva- 
tion of nebular abundances that are independent of the sys- 
tematic uncertainties of the strong line methods described 
below. In our NGC 4625 sample three H n regions were found 
to have a reliable [O ill] A4363 detection. With routines in 
the iraf nebular package we obtained T e and subsequently 
the O/H and N/O ratios, assuming that N/O = N+/0+. 
These quantities are summarised in Table [3] 

3.2 Strong Line Abundances 

Without detection of the temperature-sensitive [O ill] A4363 
line we must resort to strong line methods to determine the 
oxygen abundances. Several such methods, using a variety 
of line ratios, are present in the literature, with calibrations 
that are provided either empirically or by theoretical 
considerations. Unfortunately, it is well established that the 
various methods provide differing solutions for the chemical 
composition of star forming galaxies and Hll regions, 
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Table 3. Auroral line-derived abundances. 



ID 


[Ora] A4363/H/3 


Te (K) 


12 + log(0/H) 


log(N/0) 


(1) 


(2) 


(3) 


(4) 


(5) 


104.204 


2.0 ±1.5 


14000 ± 2700 


7.95 ±0.32 


-1.09 ±0.34 


108 


8.2 ± 1.4 


15300 ± 1100 


8.01 ± 0.09 


-1.27 ± 0.12 


205 


7.7 ±4.2 


15100 ±2500 


7.93 ±0.25 


-1.35 ±0.31 



resulting in systematic shifts of up to ~0.7 dex ( |Kewley| 
fc Ellison| [20081 |Bresolin et~al]|2009b|. By providing the 



reader with results from different strong line indicators 
we hope to detect any robust trends in the abundance 
gradient of NGC 4625, and at the same time highlight 
the limitations of our findings. Although these indicators 
were originally developed for the study of the chemical 
abundances of luminous extragalactic Hn regions, tests 
on Galactic nebulae ionised by single or only a few stars 
have excluded the presence of systematic effects ( |Kenni-| 
cutt et al.||2000| |Oey fc Shields||2000[ ). Further empirical 
tests were carried out by our team when analysing a set 
of faint outer disc H II regions in M83 ( jBresolin et al.|2009a[ ). 

We have considered the following strong line abundance 
indicators in the remainder of our work: 

1. R23 = ([On]A3727 + [Om] AA4959,5007)/H/3. This 
is one of the most widely used strong line abundance 
indicators, however it does have some limitations (e.g. 



Perez-Montero & Diaz 2005). In particular, the R23 indi- 



cator has a double valued nature, peaking with a value 
log(R 2 3)~l around 12 + log(Q/H) = 8.5 in the theoret- 



ical calibration provided by McGaugh (19911. Smaller 



R23 values correspond to both lower and higher O/H 
abundances. The turnover occurs at lower abundances, 
12 ± log(0/H) ~ 8.0, when using O/H empirical determi- 
nations from [Om]A4363 detections ( |Bresolin|[2007l [2008] ). 
Because of this degeneracy, it is essential to establish 
whether the H 11 regions under scrutiny fall into the lower or 
upper R23 branch. This is commonly accomplished with the 
help of line ratios that vary monotonically with metallicity, 
such as [Nil] /Ha and [Nil] /[On]. 

2. N2 = [Nil ] A6583 /Hq, as empirically calibrated by 
Pettini fc Pagel| fl2004| ). N2 is a monotonic function of the 
oxygen abundance and thus avoids the degeneracy problem 



associated with R23. However, Kewley & Dopita (2002 



= KD02) showed that N2 is sensitive to the ionisation 
parameter. 

3. N202 = [Nil] A6583/[Oii] A3727, a useful indicator be- 



cause it is highly insensitive to the ionisation parameter ( Do- 
pita et al.|2000[ ). We have used two separate calibrations in 



this study, the theoretical one by KD02, based on grids of 
photoionisation models, and the empirical one by |Bresolin] 
(2007 =B07), based on a sample of Hn regions with avail- 



able T e abundances. 

We note that N2 and, to a lesser extent, R23 are more 
robust against uncertainties in the atmospheric extinction 
correction compared to the N202 indicator, which is formed 
by emission lines spread across a very wide spectral range. 
N2, in particular, is also virtually insensitive to uncertain- 



ties in the interstellar reddening correction and to the po- 
tential effects of variable seeing, which could be significant 
for flux ratios that involve emission lines measured from sep- 
arately obtained blue and red spectra, as is often the case 
for [Nil] /[Oil]. 



3.2.1 Results 

In Fig. [3] we present the abundances we derived from the 
strong line indicators, plotted as a function of galactocen- 
tric distance. The three targets with T e -based abundances 
appear as gold stars. We show with black dot-dashed lines 
the weighted least squares fits to the data for both the in- 
ner (R<R 2 5) and outer (R>R 2 s) disc. We also include 
as a comparison the sample of M83 Hn regions studied 
by |Bresolin et al.| pOTOa] green dots), further discussed in 
§3.2.2| 

The triangles in Fig. [3] represent measurements of H 11 
regions in NGC 4625 from PINGS (PPAK Integral Field 
Spectroscopy Nearby Galaxy Survey, see Rosalcs- Ortega| 
|et al.|[2"010| ), which obtains spectra from the PMAS inte- 
gral field unit at the Calar Alto 3.5 m telescope. To com- 
pare the PINGS results to our own we followed two different 
approaches, either by averaging spectra from all fibers con- 
tained within equally spaced (4" separation) radial bins (red 
triangles in Fig.[3|, or by using the sum of individual fibers 
covering the spatial extent of individual Hn regions (grey 
triangles). The agreement between the abundances derived 
from PINGS and our data is excellent, as shown by the close 
match between the regressions to the average radial bins 
(red lines) and to our inner disc results (black). This can be 
quantitatively appreciated by looking at the coefficients of 
the weighted linear fits to the the radial trends of metallicity 
that are summarized in Table [4] 



Following is a brief discussion of the results obtained 
from the individual abundance diagnostics: 

-R23 - The interpretation of the R23-based abundances (top 
right-hand panel of Fig. [3| requires some important con- 
siderations. All the inner disc Hn regions in NGC 4625 
have log ([Nil] /Ha) > -0.8, which, according to the KD02 
photoionisation models, corresponds to the upper branch of 
R23. The result could be more ambiguous for the outer disc 
Hn regions, which are closer to the turnover region in the 
R23 vs. metallicity relationship, even though virtually all 
lie in the upper-branch regime, according to the |Kewley fc| 
Ellison| pOOli) criterion, log ([N 11] /Ha) > -1.1. We reach 



the same conclusion from the log ([Nil]/ [0 11]) values (rang- 
ing between —1.1 and —0.7). Therefore, we used the upper 
branch solution for the inner disc, whilst for the outer disc 
we used both the lower branch (light blue) and the upper 



branch (dark blue) calibrations of McGaugh (1991 in the 
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Figure 3. The radial oxygen abundance gradient in NGC 4625 (coloured circles), determined using various strong line methods: (a) N2, 
calibrated by |Pettini fc Pagel| p004} ; (b) R23, with upper branch calibration shown in blue, and lower branch calibration in light blue; 
(c) N202, theoretical calibration by Kewley & Dopita (2002]); (d) N202, empirical calibration byjBrcsolin ( 2007]). T e -based abundances 
have been added to all four panels as gold star symbols. We include the results obtained from measurements of the inner disc from the 
PINGS survey, including individual Hn regions studied in our sample (grey triangles) and radially averaged abundances (red triangles). 
For comparison, the M83 data points from |Bresolin et al.| l |2009a| are shown by the green dots. The lines represent weighted least square 
fits to our Subaru data for the inner and outer disc (black) and to the radially averaged PINGS data (red). 



analytical form given by Kobulnicky et al. 1999) for com- 
parison. 

Taking at face value the information obtained from the 
[Nil] /Ha and [Nil] /[On] line ratios we would be led to 
favour the upper branch solution for the whole sample of 
Hn regions, leading, in particular, to 12 + log(0/H) ~ 8.6 
in the outer disc. However, this is in evident contrast with 
the direct abundances that we measured for the three outer 
disc targets presented in Table [3] Their oxygen abundance, 
12 + log(0/H) ~ 8.0, agrees very well with the lower branch 
solution shown in Fig. [3] This inconsistency is the result of 
the still unresolved discrepancy between theoretical and di- 



rect calibrations of strong-line indicators. We note that our 
measurements for the direct abundances essentially fall in 
the transition range 12 + log(0/H) =8.0 - 8.25 between up- 
per and lower branch of the R23 calibration based on the 



P-parameter of Pilyugin & Thuan (20051, which is directly 



tied to the empirical, [O ill] A4363-based oxygen abundances. 
Consequently, we cannot obtain a reliable estimate of the 
outer disc abundances from this method, although clearly 
they would be systematically lower than those from the the- 



oretical McGaugh (19911 



N2 - We point out that this diagnostic, when applied to the 
three targets with a A4363 detection, provides abundances 
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Table 4. Coefficients of least-squares fitting to the relation 12 + log(0/H) = a — b (R/R2s). 





Abundance 


Inner Disc 


Outer Disc 




Indicator 


a 


b 


a 


b 


Our Sample 


N2 


8.72 ±0.03 


-0.17 ± 0.08 


8.38 ± 0.03 


-0.05 ± 0.02 




R23 (upper branch) 


9.02 ± 0.04 


-0.39 ±0.13 


8.38 ± 0.08 


±0.14 ±0.05 




R23 (lower branch) 






8.09 ±0.07 


-0.03 ±0.03 




N202 (KD02) 


9.08 ±0.02 


-0.25 ±0.08 


8.79 ± 0.06 


-0.14 ±0.04 




N202 (B07) 


8.63 ±0.02 


-0.27 ± 0.08 


8.28 ±0.04 


-0.10 ±0.03 


PINGS 


N2 


8.69 ±0.01 


-0.17 ± 0.03 






(Radial data) 


R23 (upper branch) 


8.94 ±0.04 


-0.28 ±0.08 








N202 (KD02) 


9.09 ±0.02 


-0.26 ±0.04 








N202 (B07) 


8.63 ±0.02 


-0.27 ± 0.05 







that are 0.3 dex higher relative to the direct method (top 
left-hand panel of Fig . [3]). This is not fully expected, since 
the N2 calibration by Pettini & Pagel (20041 relies mostly 



are also smaller (la) for N2 (-0.05 dex/R 25 ) than for N202 
(—0.10 to —0.14 dex/R,25). This agrees with the conclusion 



on extragalactic Hn regions where the oxygen abundances 
had been derived from auroral line detections However, a 



similar effect was observed in M83 by Bresolin et al. ( 2009a) 



We can speculate on the possible causes of this systematic 
difference, pointing out that the the ionization properties 
of the calibration sample can differ significantly from those 
of the faint nebulae found in the extended discs, and that 
the N2 diagnostic is sensitive to variations in the ionization 
parameter (KD02). 

N202 - The lower right-hand panel of Fig. [3] shows the 
strong line abundances calculated with the calibration of 
the [Nil] /[On] vs. O/H relation by lBresolin] ( |2007[ ). Unsur- 
prisingly, this indicator provides a good match, within the 
uncertainties, to the direct T e -based abundances, since the 
calibration was empirically derived from extragalactic H 11 
regions with auroral line detections. The well-known system- 
atic discrepancies associated with the use of different strong 
line indicators and their calibration are evident in 



The theoretical N202 calibration of Kewley & Dopita 



yields abundances that are roughly 0.45 dex higher than the 



empirical calibration of Bresolin (2007). Similar systematic 



shifts are present between N2 and N202. 



3.2.2 What do the Strong Line Abundances Reveal? 

Despite the systematic offsets between different abundance 
diagnostics noted above, the results summarised in Fig. [3] 
appear to be qualitatively consistent with each other, in 
particular showing a shallower gradient in the outer disc, 
as well as an abundance discontinuity between the inner 
and outer disc, occurring near R25. Extrapolation of the 
inner and outer disc abundance gradients suggests a discon- 
tinuity of approximately 0.15-0.2 dex, based on the N2 and 
N202 abundance indicators. As Table [4] shows, the slopes 
of the gradients depend somewhat on the choice of abun- 
dance indicator. In the inner disc, N2 produces the flattest 
gradient (— 0.17dex/R2s), about half of the value obtained 
from R23 (— 0.39dex/R25). The scatter of the data points 
is, however, quite large in the R23 case. If we look at the 
radially averaged PINGS data we find a good agreement be- 
tween R 23 and N202 (~ -0.27dex/R 25 ), still slightly larger 
(2cr) than the slope provided by N2. The outer disc slopes 



by Bresolin et al. ( 2009b I in their detailed study of NGC 300 



that abundance gradient slope values have a dependency on 
the method adopted to measure them, contrary to other 
statements in the literature, that changing the strong-line 
abundance calibration only affects the zero points. 

We include in Fig. [3] the data on the inner and outer 
H 11 regions in M83 analyzed by Bresolin et al. ( 2009a green 



dots). This allows us to compare the results we obtained for 
NGC 4625 to the only other extended disc galaxy studied in 
detail until now. As the figure shows, there is a very good 
qualitative match between the two galaxies. The change in 
the galactocentric abundance slope between inner and outer 
disc is comparable, as is the discontinuity in O/H detected 
around the isophotal radius, as well as the slope of the 
gradients themselves. The outer discs in the two galaxies 
have remarkably similar oxygen abundances, regardless of 
the metallicity indicator used for the comparison. The in- 
ner disc abundances in M83 are ~0.15 dex higher than in 
the case of NGC 4625, which is, at least qualitatively, ex- 
pected from the higher absolute luminosity of the former 
(Mb = -20.8 vs. M B = -17.2 adopting distances of 4.5 
and 9.5 Mpc, respectively). Panel (b) of Fig. [3] also suggests 
that the apparent positive abundance gradient obtained in 
the outer disc of NGC 4625 from the use of the upper branch 
of R23 could be due to small number statistics at large ra- 
dius, since the scatter of the points is similar for the two 
galaxies, and M83 does not show a positive gradient in its 
outer parts. This would imply that the R23-inferred change 
of slope between inner and outer disc of NGC 4625 would 
be made roughly consistent with what is observed with the 
other diagnostics. 



4 ANALYSIS OF THE H II REGION SAMPLE 

Earlier on (Fig. [5| we showed the familiar BPT diagnostic 
diagrams for our sample in NGC 4625, concluding that all 
objects were indeed photoionized nebulae. However, these 
plots also revealed a subset of outer disc H 11 regions which 
apparently lie off the main ionisation sequence defined by 
extragalactic giant Hn regions in the [Nil] vs. [Olll] plane, 
whilst apparently conforming to the less robust sequence 
formed in the [S 11] vs. [O ill] plane. 
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Figure 4. Properties of outer disc Hn regions, for objects with 
[Om]/H/3 > 0.25 (blue squares) and [Oni]/H/8 < 0.25 (red tri- 
angles) marked separately. Top-left: [Nll]/Ha vs. galactocentric 
distance. Top-right: [Nll]/Ha as a function of ionising luminos- 
ity. Bottom-left: ionising luminosity as a function of galactocen- 
tric distance. Bottom-right: spatial distribution of individual Hn 
regions. Inner disc Hll regions are indicated by the green dots. 



4.1 Properties of Hll Regions in the Outer Disc 

We have arbitrarily subdivided the outer disc Hll regions 
in NGC 4625 by making a cut at [Om]/H/3 = 0.25, rep- 
resented in Fig. [2] with a horizontal dashed line. The outer 
Hn regions which lie on the tight, well defined ionisation 
sequence shown in the left panel of Fig. [2] fall above our cut 
and are highlighted with central cyan dots. Those which lie 
off the ionisation sequence and below our [O m]/H/3 cut were 
marked with red dots. It is not obvious whether the latter 
deviate as a result of either a smaller [N 11] or a smaller [O in] 
emission (or both) relative to the sequence defined by bright 
Hn regions. We note that a few of the outer disc objects in 
M83, as well as a small number of nebulae in other galaxies, 
also lie below and to the left of the ionisation sequence. 

It is worth testing whether these apparently deviant H 11 
regions share any common attributes, such as ionising lumi- 
nosity, spatial distribution, metallicity and so forth. This 
is done in Fig. [4] The top-left panel shows the [Nil] /Ha 
line ratio, i.e. the N2 chemical abundance indicator, as a 
function of galactocentric radius. No differential behaviour 
between the two groups of outer disc H 11 regions (blue and 
red symbols) is evident. The two groups also appear to share 
similar ionising luminosities (upper right and lower left pan- 
els) and to be evenly distributed throughout the XUV disc 
of NGC 4625, without the presence of obvious separate sub- 
structures (lower right). There are slight indications that the 
deviant Hll regions have a systematically higher [Nil] /Ha 
line ratio than the other outer disc Hll regions, but this 
could easily be explained by the variance of the two distri- 
butions. 
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Figure 5. Histograms of the H ionising photon flux, Qo, of Hll 
regions in our sample, uncorrected for slit losses. The green and 
orange shaded areas represent Hll regions in the outer and in- 
ner disc of NGC 4625, respectively. We have marked the median 
values for the two populations (logQo =48.52 and 49.45, respec- 
tively). The blue lines show the histogram of the M83 sample 
studied by Brcsolin ct al. (2009a]l. The triangles at the top rep- 
resent reference Qo values for single hot stars (03, 06 and 09), 
from |Martins etH] | |2005[ l . 



4.2 Systematics in the Diagnostic Diagrams 

In the following we consider what might affect the line fluxes 
in such a way as to produce the observed deviation of a small 
number of outer disc Hll regions from the main ionisation 
sequence in the BPT diagnostic diagram. 



4-2.1 Stochastic Variations 

The direct observation of the ionising population of nebu- 
lae in outer galaxy discs are extremely difficult, thus the 
properties of the associated stellar clusters have so far been 
obtained from studies of their Ha and UV luminosities ( |Gil| 
de Paz et al] [20051 |2008| jGoddard et al.|[20T0| , together 



with their IR emission ( Dong et al.|2008 l. The total ionising 



photon flux of individual clusters can be easily estimated 
from the Ha luminosity, using the conversion given by |Ken-| 
< [l998l >. In Fig. [5] we show histograms of the hydrogen 



nicutt 



ionising photon flux for both the inner (orange) and outer 
(green) disc Hll regions, along with their median values. 
As a reference, we also show the hydrogen ionising fluxes 
of 03, 06 and 09 dwarf stars, taken from Martins et al. 
( |2005[ ). The luminosities that we measure do not account 
for slit losses, which are expected to be smaller for the outer 
disc H 11 regions, given that a greater percentage of their flux 
will fall through the 1" slits compared to the brighter and 
larger inner disc Hll regions. We estimated the slit losses by 
measuring H 11 region fluxes on an Ha image of NGC 4625 
provided by the SINGS survey (Kennicutt et al. 20031. In 



the case of the faintest Hll regions losses are in the 5-10% 
range, while for the brighter inner H 11 regions they increase 
to 20-70%. 

Figure [5] confirms what other authors already found, 
namely that that the inner disc sample of Hll regions is 
brighter than the outer disc population by at least one or- 
der of magnitude (more if slit losses were considered). The 
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ionising flux of typical outer disc H II regions corresponds to 
that of a single late-type O sta r (see|Gil de Paz et al.|2005 



Thilk er et al.|2005| |Dong et al.|2008[|Goddard et al.|20To [ 

Consequently, it is necessary to account for the stochastic 
nature of the upper Initial Mass Function (IMF) of the ion- 
izing clusters. An ionising population drawn from a stochas- 
tically sampled IMF, as opposed to a fully sampled one, will 
tend to have a lower effective temperature, due to the pref- 
erential depletion of the rarest, most massive main sequence 
stars. This has important consequences for the total ionizing 
flux output, as well as for the hardness of the ionizing radia- 
tion. In particular, as a result of softer spectral energy distri- 
butions, line ratios involving both high- and low-ionization 
potential species could be significantly affected. 

4.2.2 Age effects 

As the ionizing population of an H 11 region ages the most 
massive stars die first, effectively truncating the mass func- 
tion of the remaining stellar population. Not only could this 
alter the emission line fluxes, but we can also think of this 
effect as a proxy for stochastic variations which, as we have 
just seen, also act in the send of removing the high mass end 
of the stellar population. 

We make a simple attempt to model such age effects by 
using the itera program ( Groves & Allen 2010 1, along with 



the starburst models of Levesque et al. (20101, which adopt 



a Salpeter IMF and the Starburst99 (Leitherer et al.||1999| ) 
stellar spectral energy distributions with standard mass loss 
rates. In Fig. [6] we show a set of models on the [Om]/H/3 
vs. [Nil] /Ha diagnostic diagram, calculated for a range of 
metallicities (Z = 0.2, 0.4, and 1.0 Zq), ionisation parame- 
ters (Qo = 2 x 10 7 , 8 x 10 7 , and 4 x 10 8 ) and ages (t = 0, 3, 5, 
6 Myr) . The H 11 region samples introduced in ij2] are repre- 
sented by the grey dots, whilst our targets in NGC 4625 
are shown either with cyan ([Oiii]/H/3 > 0.25) or green 
([O m] /H/3 < 0.25) symbols. 

The zero-age models in Fig.[6]lie on top of the ionisation 
sequence defined by the main extragalactic sample of giant 
H 11 regions. The models show that, as the H 11 regions evolve 
with time after the initial burst of star formation, the [O ill] 
emission rapidly decreases, and the good match with the 
observational data quickly disappears (see earlier results by 
Bresolin et al.|1999 ). The models also depart from the main 
ionisation sequence to lower [Nil] /Ha values as the metallic- 
ity decreases. The best fit to our deviant Hn regions (green 
dots) is obtained with the 0.4 Zq models and ages between 
3 and 5 Myr, depending on the ionisation parameter. 

The model comparison offers some clues as to why some 
Hn regions appear to deviate from the main ionisation se- 
quence. Our interpretation relies on the fact that, as men- 
tioned earlier, the evolution of the ionising output of low- 
luminosity Hn regions with age can mimic the effects due 
to stochastic variations in the number of ionizing stars. One 
might ask why so few H 11 regions have been seen with simi- 
lar properties in previous studies. The reasons are two-fold. 
Firstly, most of the previous studies have focused on inner 
disc H 11 regions, which tend to have higher metallicities. Ac- 
cording to our models, these objects would be nearly super- 
imposed onto the main ionisation sequence. Secondly, there 
is a selection effect at work, since in the inner discs the 
brightest and youngest Hn regions are preferentially stud- 




1.0 -0.5 

log([NII]H6548,6583/Ha) 

Figure 6. Models in the [Oni]/H/3 vs. [Nil] /Ha diagnostic dia- 
gram. The grey dots represent the same extragalactic H II region 
samples included in previous figures. Cyan symbols are used for 
outer disc Hll regions in NGC 4625 for which [Om]/H/3 > 0.25, 
whilst green dots are used for those with [OlIl]/H/3 < 0.25. The 
red line represents the upper boundary for photoioniscd nebulae 
defined by Kcwlcy et al. ( 2006 ). Black lines link models with iden- 
tical ionisation parameter, whilst the dotted lines link models of 
the same age. 



ied, also as a result of source confusion and crowding near 
the central regions and along the spiral arms, while in the 
outer disc we attempt to sample as many Hll regions as 
possible, even the faintest ones, thus probing to older ages. 
This conclusion could be tested with spectroscopic observa- 
tions of larger samples of inner disc Hll regions, reaching to 
fainter emission levels. 

We estimated the ages of the Hll regions in our sam- 
ple by looking at the UV and Ha fluxes measured by God- 
dard et al.| ( |2010[ ). In Fig. [7] we plot the FUV/Ha ratio as 
a function of the FUV luminosity. We have marked with 
filled circles those objects for which an Hll region matches 
the position of an identified FUV object. For the remaining 
cases (open circles) we have measured the total Ha emission 
coming from a region that matches in size the UV-identified 
objects. We include Starburst99 models calculated for a sin- 
gle burst of star formation at three different cluster masses 
(10 2 , 10 3 and 10 4 Mq). Dashed lines connect models of iden- 
tical age. It can be seen that the cluster ages predicted by 
the models are roughly consistent with the ages of the it- 
era models at which deviations in the BPT diagram occur. 
However, the data are unable to break the age-stochasticity 
degeneracy, and thus we cannot differentiate between the 
two effects. 

Fig. [7] also provides estimates for the embedded cluster 
masses, indicating values ranging between 10 3 and 10 4 M . 
Such low cluster masses are in the regime that is under the 
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Figure 7. FUV/Ho ratio as a function of FUV luminosity for 
XUV disc objects in NCG 4625 as measured by Goddard ct al. 
( |2010| . Objects for which an H II region in our sample falls under 
the area of a UV-detected object are shown as filled circles. Open 
circles are used in cases where we have instead measured the total 
Ha emission coming from a region that matches in size the UV- 
identificd objects. Grey lines show theoretical Starburst99 models 
for three cluster masses (10 2 , 10 3 and 10 4 Mq). Dashed lines 
connect models of constant age. 



influence of stochastic variations in the upper end of the 
stellar mass function, as already pointed out by |Gil de Paz| 
et al.l (12005b. 



4.2.3 Geometry 

In an Hn region with a small number of ionising sources, 
such as those which have a stochastically selected popula- 
tion of massive stars, the spatial distribution of the ionis- 
ing sources relative to the nebular gas can affect the ioniz- 
ing properties of the cluster. |Ercolano et al.| ( |2007| ) showed 
that geometric effects could produce theoretical H n regions 
which lie off the main ionising sequence in the BPT diag- 
nostic diagrams, as well as affect the nebular abundance 
determination, especially at low metallicities. However, no 
particular model from this paper was found to reproduce 
our observations. 



4.2.4 Anomalous N/O ratios 

In the outer discs of galaxies the star formation rate is par- 
ticularly low, typically only a few percent of that in the inner 



discs, moreover spread over a much larger area ( Gil de Paz 
|et al.||2005| |Goddard et al.||2010| ). In proposing that small 
stellar clusters are unable to form massive stars, |Pflamm^| 
Altenburg & Kroupa ( 2008 1 suggested that this might ex- 



stars in the outer discs, and therefore of core-collapse su- 
pernova explosions, we would expect a reduced enrichment 
of oxygen relative to nitrogen, since oxygen is produced by 
the most massive (M > 8 Mq) stars only, while nitrogen 
is also produced by intermediate-mass stars. Could such an 
effect be responsible, at least in part, for the anomaly that 
we found in the BPT diagram? 

The left-hand panel of Fig. [8] shows that the deviant 
H II regions in the BPT diagnostic diagram (green dots) are 
shifted to the left relative to the sequence formed by bright 
Hll regions also in the R23 vs. [Nii]/[Oii] plane. The shift 
appears to occur at constant [Nii]/[Oii] ratio, at variance 
with what one would expect if the gas was lacking in oxygen 
relative to nitrogen, which would move the objects towards 
the top of the diagram. 

4.3 Reliability of the Strong Line Abundances 

We still need to consider the fact that the systematic offsets 
found in a subset of the outer disc H 11 regions in NGC 4625 
could, in principle, affect the abundances determinations 
based on strong line indicators. Would these indicators, de- 
veloped for bright extragalactic ionized nebulae, still be ap- 
plicable to much fainter Hll regions, where the population 
of massive stars is stochastically determined? 

In Fig. [8] we show model predictions and empirical data 
for R23 as a function of the [Nil] /[On] ratio (left panel) 
and of the [Nil] /[On] ratio (right panel). In the latter case 
we plotted lines of constant metallicity, adopting both the 
upper (black lines) and lower (grey) branch of the theoretical 
calibration by |McGaugh| ( |1991[ ). 

The left-hand panel of Fig. [8] gives us two major in- 
sights. Firstly, it confirms our earlier finding (i 4.2.2 1 that 



plain the edge of Ha emission in galactic discs and the lack of 
massive stars beyond this radius. With a dearth of massive 



the line emission properties of ageing Hll regions appear 
to replicate those of the Hll regions in our sample which 
were found to lie off the main ionisation sequence in the 
BPT diagnostic diagram (green dots). Secondly, it shows 
how an ageing population of Hll regions might affect the 
determination of the chemical abundances. The itera mod- 
els suggest that as an Hll region ages, its R23 value rapidly 
decreases, resulting in a shift almost parallel to the lines of 
constant metallicity (determined using the N202 diagnos- 
tic). The [Nil] /[On] ratio also changes, simulating a varia- 
tion in the inferred metallicity. Therefore, the models give 
us an indication of how reliable the metallicity indicators 
become for an ageing Hll region population. A spread of 
ages of, say, between 3 and 6 Myr, as suggested by Fig.[7]for 
NGC 4625, could result in an uncertainty in the abundance 
measurement of roughly ~ 0.15 dex, based on the N202 
diagnostic. 

The right-hand panel of Fig. [8] helps us to identify how 
robust the R23 indicator may be against such an ageing H 11 
region population. For a metallicity of 0.4 Zq the Hll re- 
gions appear to evolve almost parallel to the lines of constant 
metallicity, although there may be some divergence with the 
lower branch calibration for the oldest models. However, for 
the 1 Zq metallicity case the H 11 region evolution appears 
to be problematic for the models, since large changes in the 
measured abundance, of the order ~0.25 dex, appear to oc- 
cur, together with discrepancies for models of the same age 
but lower ionisation parameters. This suggests that for an 
ageing population of H 11 regions with a spread of ionisation 
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Figure 8. Left: log R23 vs. log [Nn]/[Oll]. We show lines of constant metallicity based on the N202 calibrations of |Kewley fc Dopita| 
(2002, black dashes) and|Brcsolin (2007] grey dashes). Right: log R23 vs. log [Om]/[Olrj. We show lines of constant metallicity, based on 
the upper (black lines) and lower (grey) branches of the R23 calibration of McGaugh (1991). In both panels the grey dots show the same 
sample of H II regions plotted in previous figures. Cyan dots represent H II regions from our NGC 4625 sample for which [O IIl]/H/3 > 0.25, 
whilst green dots represent those with [Olll]/H/3 < 0.25. We have also plotted ITERA nebular models using the same symbols as in Fig.RH 



parameters the R23 indicator could represent an inaccurate 
measure of the oxygen abundance. 



5 DISCUSSION 



5.1 Breaks in Abundance Gradients 



In §3.2| we showed that the strong line abundance indica- 
tors used in this study gave qualitatively similar results. All 
showed a discontinuity in the abundance gradient at a galac- 
tocentric distance between 0.7 and 1.1 R25, and a shallower 
gradient in the outer disc. 

While the abundance solution provided by the use of the 
R23 indicator remains somewhat dubious, due to the uncer- 
tainty in the choice between upper and lower branch cali- 
bration, the N2 and N202 abundance diagnostics are mono- 
tonic with abundance and are unlikely to produce false dis- 
continuities or gradient changes. In addition, Brcsolin ct al. 
| |2009a| ) tested the reliability of N202 to measure relative 
abundances in both the inner and outer disc of M83. 



In both NGC 4625 (this work) and in M83 ( |Bresolin 
et al.|2009a l we found a discontinuity in the abundance gra- 



dient, with a shallower slope beyond R25. Goddard et al.| 
(20101 also found similar radial UV and Ha profiles in the 



two galaxies. The FUV profile is smooth over the transition 
between inner and outer disc, becoming shallower beyond 
R25. Conversely, the Ha profiles show a sharp truncation at 
R25, in correspondence with the edge of the star forming 
disc. The origin of these effects lies possibly in the differing 
star formation properties between inner and outer parts of 



the discs we have analyzed. The extremely low Ha emission 
in the extended discs indicates a very low star formation rate 
(Gil de Paz et al.||2005| iGoddard et al.||2010h. The lack of 



massive stars results in a smaller chemical enrichment from 
supernova events, which would eventually lead to the oxy- 
gen abundance discontinuity observed between the inner and 
outer discs. G07 have argued that both a low-level, but con- 
tinuous, star formation activity and an episodic one could 
be invoked to explain the fact that the oxygen abundances 
are low, but not extremely so, in the extended discs of spi- 
rals (see also Bresolin et al.|[2009a ). Bresolin et al. (2009a I 
have remarked on the similarity between XUV discs and low- 
surface brightness galaxies in terms of chemical abundances 
and enrichment. They also pointed out that gravitational in- 
teractions with companion galaxies could generate, via gas 
stripping, peculiarities in the abundance gradients, such as 
those observed in the outer disc of M83, which has likely 
interacted with NGC 5253 around 1 Gyr ago. The situation 
appears to be similar in the case of NGC 4625, which has 
been in interaction with its companion NGC 4618 and, pos- 
sibly, a fainter galaxy ( Gil de Paz et al.|20 05). Perhaps these 
interactions play an important role both in the activation of 
the star formation in the extended discs, and in the genera- 
tion of the abundance breaks and/or flattened gradients. 

5.2 HII Region populations in XUV discs 

We found a significant number of H 11 regions in NGC 4625 
that appear to depart from the main ionisation sequence es- 
tablished by studies of bright, inner disc ionized nebulae. A 
similar behaviour is detected (Fig. [2| for a small fraction of 
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H n regions in other galaxies. The most likely reasons for this 
anomaly have been identified in stochasticity or the effects 
of an ageing H II region population. Since the outer disc H II 
regions contain at most only a few O-type stars, the source of 



ionizing flux is subject to stochastic variations (Gil de Paz 


et al. 2005 Thilker et al. 


2005 Goddard et al. 20101. By 


using measurements from 


Goddard et al. (20101, in combi- 



nation with published photoionisation models, we also found 
that an ageing Hn region population might equally be re- 
sponsible for the spread we observe in the BPT diagnostic 
diagram. These two effects are indistinguishable, since an 
evolving Hn region loses first the most massive stars, which 
are also those that are most affected by stochastic variations. 

We examined the possibility that geometry or a non- 
standard N/O ratio might explain the emission line proper- 
ties of the deviant Hn regions. However, we concluded that 
these effects are unlikely to produce the observed trends, 
and both age and stochastic arguments appear to be more 
robust. It was never our intent to produce a full description 
of these effects, but only to highlight the fact that outer 
disc Hn regions are significantly different from their inner 
disc counterparts, and this difference has not been fully in- 
vestigated yet. In the future it will be important to investi- 
gate the presence of H II regions in other galaxies with emis- 
sion line properties similar to what we found in NGC 4625, 
and carry out new comparisons to models generated with a 
stochastically populated IMF. 



6 CONCLUSIONS 

In this study we have analysed the spectra of 34 H II regions 
at a range of galactocentric radii in the XUV disc galaxy 
NGC 4625. Below we summarise our main findings: 

1. Strong line abundance indicators consistently show a 
shallower abundance gradient in the outer disc compared 
to the inner disc. This is accompanied by a discontinuity 
in the abundance gradient that occurs in proximity of the 
optical edge of the galaxy. These trends are very similar to 
those seen in M83. 

2. The three T e -based abundances in the outer disc predict 
a metallicity of 12 + log(0/H) ~ 7.96 ± gjf (~ 0.2 Z@) in 
the outer disc. 

3. We identified a subset of outer disc Hn regions which 
do not conform to the ionisation sequence defined by the 
bright inner disc Hn regions in the [Oiii]/H/3 vs. [Nil] /Ha 
diagnostic diagram. These objects do not appear to display 
any unique properties in terms of position within the XUV 
disc or ionising flux. 

4. Stochastic variations in the high- mass end of the embed- 
ded stellar population or an ageing Hn region population 
seem to be the most likely effects capable of producing the 
observed emission line properties. However, it has not been 
possible to distinguish between these two effects. 

5. The strong line abundance indicators used in this paper 
appear to be less reliable if the Hn region population is 
ageing, even over the short 10 Myr timescale of the ionising 
stars. This effect could be as much as ~ 0.15 dex for the 
N202 indicator, and even more for R23. 

We hope future studies of outer disc Hn regions will 
be able to further probe this regime of stochastically de- 



termined ionisation sources. It is only through the contin- 
ued study of the abundance gradients in a variety of XUV 
discs that we may hope to better understand the connec- 
tion between star formation and chemical enrichment of the 
interstellar medium in these low gas density environments. 
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